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THERMAL STRUCTURES: FOUR DECADES OF PROGRESS 


Earl A. Thornton, Professor 
Department of Mechanical and Aerospace Engineering 
University of Virginia 
Charlottesville, Virginia 22901 


Abstract 

Since the first supersonic flight in October 1947, the 
United States has designed, developed and flown flight 
vehicles within increasingly severe aerothermal 
environments. Over this period, major advances in 
engineering capabilities have occurred that will enable the 
design of thermal structures for high speed flight vehicles 
in tne twenty-first century. This paper surveys progress 
in thermal-structures for the last four decades to provide a 
historical perspective for future efforts. 

i n tr od uctio n 

The design of structures for winged flight vehicles 
that fly through the earth's atmosphere, either to and from 
space or in sustained flight, poses severe challenges to 
structural designers. Major components of the challenge 
are to select materials and design structures that can 
withstand the aerotheimal loads of high speed flight. 
Aerothermal loads exerted on the external surfaces of the 
flight vehicle consist of pressure, skin friction (shearing 
stress), and aerodynamic healing theat flux). Pressure and 
skin friction have important roles in aerodynamic lift and 
drag, but aerodynamic heating is the predominant 
structural load. Aerodynamic heating is extremely 
important because induced elevated temperatures can 
affect the structural behavior in several detrimental ways. 
First of all, elevated temperatures degrade a material's 
ability to withstand loads because elastic properties such 
as Young's modulus are significantly reduced. Moreover, 
allowable stresses are reduced, and time-dependent 
material behavior such as creep come into play. In 
addition, of course, thermal stresses are introduced because 
of restrained local or global thermal expansions or 
contractions. Such stresses increase deformation, alter 
buckling loads and change aerotbermodynamic behavior. 

Today, the advent of the National Aerospace Plane 
offers structural engineers new challenges for the design of 
thermal structures for high speed flight. A theme of this 
paper is that the progress in research and development of 
thermal structures and related technologies over the last 
four decades provides the foundations to meet the new 
challenges. 

One justification for a survey paper is expressed 
very well by the opening of the "House divided" speech 
given by Abraham Lincoln in 1858, "If we could first 
know where we are, and whither we are tending we could 
better judge what to do, and how to do it". Tne speech, 
iven to a political convention in Springfield, Illinois in 
une of that year, referred to an altogether different 
subject, but surely the words apply in the present context. 
As our nation looks to high speed flight vehicles for the 
next century, we should review past efforts to provide a 
historical perspective for the future, 

Aerospace technology has advanced far in a very 
short time. These advances are particularly notable for 
high speed flight vehicles. This paper attempts to provide 
a historical perspective for thermal structures by 
describing the evolution of thermal-structures technology 
from the early 1950's to the technology of the early 199(rs. 


The early 1950's was selected as the reference point 
because the first manned supersonic flight in 1947 
stimulated a period of intense research on high 
temperature structures. We are fortunate today to have 
access to the papers, reports and books that describe 
research of that era. Unfortunately, although many 
researchers of that period are alive today, much of their 
personal experience is not available to younger researchers 
because of the "generation gap" that currently exists in the 
aerospace establishment. Thus for all of us interested in 
thermal structures, there are lessons to be learned from 
surveying the progress of the last forty years 

Evolution, of Thermal Structures 
From World War II to Sputnik 

The need to understand aerothermal loads ami the 
design of thermal structures have their origins in the late 
1940’s. In World War II. airplane speeds had become high 
enough for compressibility phenomena to have a significant 
role in performance. Transonic phenomena were not well 
understood, and over a period of years the phrase "sound 
barrier" came into use. The need for a transonic research 
airplane was recognized during the war, and in 1944 the 
design development of the Bell X-l program was initiated 
(ref.l). The X-l proved tu be enormously successful, and 
the llight of Captain Charles E Yeager on October 14. 
19^7 proved beyond doubt that manned aircrab could fly 
faster than the speed of sound. An advanced version ot the 
aircraft, the X— IB shown in Fig 1, flew several research 
missions for NACA to study aerodynamic heating effects. 
The original X-l aircraft as well as the advanced version 
used aluminum construction throughout. Measured skin 
temperatures are shown for a NACA research mission 
flown in January 1957 at Mach 1.94. Note that skin 

temperatures are low, less than 200 °F. Thereafter, 
supersonic flight speeds increased rapidly, and the need for 
considering aerodynamic beating became evident. The 
difficulties presented by high temperatures accompanying 
flight at supersonic speeds became known as the "thermal 
barrier". For about ten years, the "thermal barrier" or 
"thermal thicket" caused' concern that large structural 
weight increases would be required to keep material 
temperatures within allowable values. Subsequently 
researchers (refs. 2—3) found that these concerns did not 
materialize because the problems were overcome through 
research and development of effective thermal structures. 

After the first supersonic flight, research and 
development of high speed aircraft intensified. A contract 
for the design, development and construction of two X-2 
swept wing, supersonic research aircraft was awarded to 
the Bell Aircraft Corporation in 1947. The X-2 was the 
first aircraft structure designed for aerodynamic heating 
fref.4). Until the X-2, speed s had not been high enough 
tor the structure to be affected adversely by aerodynamic 
heating. For increased strength at elevated temperatures, 
the fuselage was constructed from K-Monel, and 
aerodynamic skins used stainless steel. A drawing shown 
in Fig. 2 illustrates the X-2 structure. The X-2 became 
the first research aircraft to explore the "thermal thicket" 
with speeds above Mach 2.5 On September 27. 1956 the 
X-2 achieved its maximum speed of Mach 3.2. but 



unfortunately the aircraft went out of control and the 
pilot, Captain Milburn G. Apt. was killed. 

During this period theoretical and experimental 
studies were in progress to develop the technologies needed 
for high speed aircraft. Considerable research was 
underway in compressible flow, thermal structures and 
materials. Excellent descriptions of these efforts are 
available through survey articles and collections of papers 
in books describing thermal-structures conferences. The 
Applied Mechanics Review article (ref. 5) of November 
1955 by Professor N. J. Hoff of Polytechnic Institute of 
Brooklyn describes high temperature effects in aircraft 
structures with emphasis on the effects of temperature on 
buckling and creep. In December, 1955 Professor R. L. 
Bisplinghoff presented the Nineteenth Wright Brothers 
Lecture to the Institute of Aeronautical Sciences in 
Washington, D. C. The paper (ref 6). including discussion 
by several prominent researchers, was published in April 
1956. Professor Bisplinghoff s paper gives a comprehensive 
review of structural considerations for high speed flight 
and an excellent account of design and analysis practices. 
A book (ref.7) published by the Advisory Group for 
Aeronautical Research and Development (AGARD) of the 
North Atlantic Treaty Organization (NATO) in 1958 
describes effects of high temperatures on aircraft structures 
caused by aerodynamic heating. The book, edited by 
Professor* Hoff, contains Hi articles written by U. S. and 
European authors. 

The next major flight program that stimulated 
thermal structural research was the X— 15. The X-15 had 
complex origins including the prewar and postwar work of 
German scientists Eugene Sanger and Irene Bredt who, in 
1944, outlined a hypersonic, rocket-propelled aircraft. 
The evolution of their ideas which contributed to the 
development of the X-15 is described in the award 
winning paper by Richard P. Hallion (ref. S). Further 
descriptions of the X-15 program aie given by Hallion in 
Ref. 9. and by NASA Langley research scientist John V. 
Becker in Ref. 10. The paper written by John Becker was 
presented in Bonn, West Germany in Dec 1968. At this 
meeting he accepted the Eugene Sanger medal awarded by 
the German Society of Aeronautics and Astronautics to 
honor the success of the X-15 program. 

On June 24. 1952 the prestigious NACA Committee 
on Aeronautics charged the agency to study problems of 
manned and unmanned flight at altitudes between 12 and 
50 miles and speeds of Mach 4 to Mach 10. By 1954 the 
NACA Langley Laboratory (now NASA's Langley 
Research Center) had formed a hypersonic study team 
with John Becker as the chairman. According to Becker 
(ref.10) the unprecedented problems of aerodynamic 
heating and high-temperature structures appeared to be so 
formidable that they were viewed as "barriers" to 
hypersonic flight. Nevertheless, the NASA group evolved 
a baseline design that closely resembled the ultimate X-15 
configuration. Three X-15’ aircraft were built by North 
American Rockwell. 

A thick-skinned, heat-sink approach was adopted 
to suit the short duration missions of the X-15. A typical 
research mission lasted 10-12 minutes (ref.ll). Surfaces 
exposed to aerodynamic heating were made of Inconel X, a 
nickel alloy. Internal structures not subject to high 
temperatures were made of titanium. Skin temperatures 

were designed for a maximum of 1200 °F. Figure 3 shows 
maximum temperatures experienced in an X-15 flight. 

On October 4. 1957 the Soviet Union orbited 
Sputnik 1, the world's first artificial satellite. 


This event changed the nation's priorities for high speed, 
high altitude flight making the X-15 program vital to 
America's national prestige. The X-15 program served 
the nation well accomplishing 199 missions between 1959 
and 1968. The X-15 was the first, and to date, the only 
manned vehicle capable of flying atmospheric missions at 
Mach 5 for altitudes of 100.000 feet or higher. It nude 
many contributions to the understanding of hypersonic 
flight including the design of thermal structures. Later 
sections will refer to some of these contributions. 

From World War II to Sputnik, the nation's 
research programs were focused on high speed flight of 
aircraft. In thermal structures, the NACA Langley 
laboratory made significant contributions. Many of the 
Langlev research efforts were led by Richard R. Heldenfels. 
After his retirement, he describes (ref. 12) research 
conducted at Langley from 1948 to 1958 on structural 
problems caused by aerodynamic heating. 

After Sputnik 

After Sputnik, the nation's high speed flight 
research program was broadened to include a major 
emphasis on manned space flight. The National 
Aeronautics and Space Act of 1958 created NASA, and 
NACA became NASA on October 1, 1958. 

The expanded scope of the research is reflected in 
the I960 paper by Heldenfels (ret 3). The paper describes 
proposed re-entry structure.^ <is well as space vehicles and 
space structures. The proceedings (ref. 13) of a conference 
held in Cambridge. Mas^u nu>etts on July 25. 1961 
contains 13 papers describing research on thermal 
structures for manned and unmanned re-entry vehicles. 
The conference was concerned primarily with thermal 
protec-: ion systems for lifting vecbicles. Two types ol 
thernul protection systems are described: 

"cool-structure" and "hot-structure" approaches, in the 
cool-structure approach, the load-bearing structure is 
insulai ed from high temperatures by an external heat 
shield. In the hot-structure approach, the load-bearing 
structure operates at nearly skin temperatures. 

A NASA— University Conference on the Science 
and Technology of Space Exploration was held in Chicago. 
Illinois on November 1-3, 1962. The conference 

proceedings describes research on structures for launch 
vehicles, winged aerospace vehicles and planetary entry 
vehicles. The paper by Mathauser (ref. 14) gives a good 
account of Langley ’ research on thermal-structural 
problems for winged vehicles. 

Types of winged vehicles under consideration bv 
NASA are shown in Fig. 4. These include a research 
airplane (the X-15), a re-entry glider, and a large 
hypersonic aircraft. The re-entry glider is representative 
of a X-20 Dyna-Soar vehicle that was to be launched with 
a ground-based booster. The hypersonic aircraft was to 
possess horizontal take-off capability and be capable of 
sustained hypersonic flight. All three vehicles utilize 
radiation-cooled hot structures. Radiation equilibrium 
temperatures (See Eq. (13)) estimated for a re-entry glider 
are shown in Fig. 5. 

The re-entry glider and hypersonic aircraft of the 
Mathauser paper are hypothetical vehicles used in NASA's 
fundamental research studies. However, the Dyna-Soar 
project sponsored by the Air Force did lead to the final 
design of the Boeing X-20. • The events leading to the 
Dyna-Soar program, its evolution, and subsequent demise 
are described by Miller (ref. 4) and Hallion (ref. 8). The 
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Boeing development program began with a contract award 
in November 1959. The X-20 was designed to provide a 
piloted, maneuverable vehicle for conducting experiments 
in the hypersonic and orbital flight regime. However, the 
X-20 program was cancelled in December 1963 before the 
first vehicle was completed. The Dyna-Soar program 
accelerated progress in several technologies that ultimately 
were applicable to the space shuttle. A review of the 
Dyna-soar winged spacecraft technology appears in the 
1961 paper by Yoler (ref. 15). The structural design 
utilized a Rene 41 nickel super alloy primary structure, a 
Columbian alloy heat shield, a graphite and zirconia nose 
cap and molybdenum alloy leading edges. 

Concurrent with these activities, the U. S.'s 
manned space flight program accelerated rapidly. The 
effort began with the Mercury program including the Alan 
B. Shepard, Jr. suborbita) flight on May 5. 1961 and 
America's first orbital flight by John H. Glenn. Jr. on 
February 20. 1962. The .effort continued with the Gemini 
program and the first two-man fight by John VV. Young 
and Virgil I. Grissom on March 23, 1965. The Gemini 
program achieved the first rendezvous and docking in 
space and the first American "space walk". The Apollo 
program began with the October 11-22. 1968 llight by 
Walter M. Shirra. Jr.. Donn F. Eisele and R. Walter 
Cunningham. The program reached its zenith with the 
historic fight of Apollo 11 by Neil A. Armstrong, Edwin 
E. Aldrin, Jr. and Michael Collins and the first lunar 
landing on July 20. 1969. The lunar program concluded 
with Apollo IT making the sixth and last lunar landing in 
December 1972. America's first Earth-orbiting space 
station. Skylab, was launched atop a Saturn V booster on 
May 14. 1973. Three, three-manned crews visited the 
space station with the last mission returning to Earth in 
February 1974. The rendezvous and docking of an Apollo 
spacecraft with a Russian Soyuz craft in Earth orbit on 
July IS. 1975 closed out rhe Apollo program. 

All of the manned spacecraft missions used 
blunt— body re-entry vehicles and ablative heat shields to 
dissipate aerodynamic heating. These ballistic or 
semi-ballistic vehicles had small lift to drag ratios 
permitting only limited maneuverability. 

The lack of a follow-on winged research fight 
vehicle to succeed the X-15 forced thermal-structures 
research in the 1960's into a period of more fundamental 
rather than applied research. The 1962 paper (ref. 14) by 
Mathauser indicates that research had begun at Langley on 
hypersonic structures. A 1966 paper by Heldenfels (ref. 16) 
discusses in more detail structural prospects for hypersonic 
vehicles. Figure 6 shows isotherms on a hypersonic vehicle 
assumed to be cruising at Mach S at 88,000 ft. 
Temperatures shown are radiation equilibrium 
temperatures expected during a typical fight of one to two 
hours duration. The paper discusses thermal-structural 
designs for fuselage tanks for liquid hydrogen, wings and 
engine structures. Hydrogen fuel-cooled structures for 
engines and passive hot structures of high temperature 
materials for airframes are described. The proceedings of a 
conference held at Langley in November, 1971 contains 
several papers describing basic hypersonic vehicle research. 
A paper by Anderson and Kelly (ref. 17) reviews the 
technology base for propulsion structures, primary- 
structures and liquid hydrogen tanks. A recent paper by- 
Shore (ref. 18) reviews research on convectively cooled 
structures in the 1960's and 1970's. 

The Space Shuttle 

The Space Shuttle resulted from a perceived need in 
the 1960‘s for a logistical spacecraft to support orbital 


space stations. However, after the lunar landing in 1969, 
NASA recognized that funds would not be available to 
support both the Shuttle and a space station. Justification 
for the Shuttle shifted from space station support to its use 
as a substitute for expendable launch svstems. With 
strong support from the Department of Defense, the 
preliminary analysis phase began in February 1969 with 
contracts to Lockheed, General Dynamics, McDonnell 
Douglas and North— American Rockwell. Rockwell 
eventually became the prime contractor, and construction 
of the actual Shuttle orbiter began in June 1974. Rockwell 
completed this vehicle, the Enterprise, in September*l976. 
Difficulties with the main engines and thermal protection 
system delayed the second shuttle, the Columbia until 
1981. Piloted by astronauts John W. Young and Robert 
L. Crippen, Columbia completed the Shuttle's first orbital 
flight on April 14, 1981. Further details of its history are 
given by Hallion (ref.8). 

The orbiter basically has a conventional 

skin— stringer alumninum aircraft structure. The design ol 
the thermal protection system (TPS) had the requirement 

of keeping structural temperatures less than 350 °F. The 
thermal protection is composed of two type> of reusable 
surface insulation (RSI) tiles. The RSI tiles covering thi j 
orbiter are made of coated si lea fiber. The two type- differ 
only in surface coating to provide protection for difieieui 
temperature environments. The low— temperature 

insulation (LRSI) consists of S-inch square silca tiles and 
covers the top of the vehicle where temperatures are less 

than 1200 °F. The high temperature insulation ( H RSI ) 
cons sts of 6— inch square tiles that cover the bottom and 
some leading edges of the oi biter where temperatures are 

below 2300 C F. Reinforced carbon— carbon (RCC) is used 
for the nose cap and wing leading edges where 

temperatures are above 2300 °F. Flexible reusable surface 
insulation (FRSI) is used at locations where temperatures 

are less than 700 °F. The distribution of the tiles is' shown 
in Fig. 7. and orbiter isotherms for a normal flight are 
shown in Fig. S. Further details of the orbiter are given in 
Ref 19. The launch of the first shuttle flight was delayed 
by problems with the TPS, and the article (ref. 20) by 
Langley researchers Paul A. Cooper and Paul F. Holloway 
explains the problems and their resolution 

The National Aerospace Plane i \-30t 

In his State of the Union address in February 19SG. 
President Reagan made the initial announcement of 
research on development of the National Aerospace Plane 
(NA3P). The November 1986 article by Robert L. 
Williams (ref. 21) describes the program objective "to 
develop, and then demonstrate in an experimental flight 
vehicle (the X-30), the requisite technologies to permit the 
nation to develop both military and civil vehicles capable 
of operating at sustained hypersonic speeds within the 
atmosphere and/or operating as space launch vehicles for 
delivering payloads into orbit" Under joint sponsorship 
by NASA and DOD, technology development began in 
19S6, and by 1988 McDonnell Douglas, General Dynamics 
and Rockwell had won separate contracts to proceed with 
development of the airframe. Pratt and Whitney and 
Rocketdyne had won contracts for propulsion system 
development. A consortium of these contractors is 
participating in a 30— month materials and structures 
augmentation project (ref. 22). 

The competitive nature of the development 
program and security issues have limited pnbliuii sons 
describing thermal structures for the X-JU. However, at 
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the 30th SDM conference, David A. Ellis of the McDonnell 
Douglas Corporation presented an overview paper (ref. 23) 
describing the concepts under consideration as well as the 
types of trade studies being conducted. Figure 9 shows 
representative surface temperatures for a typical NASP 
hypersonic condition. Candidate primary structure design 
concepts include hot structures, actively cooled structures 
and thermally protected structures. The design of the 
center fuselage is particularly challenging because the 
external skin is subjected to aerodynamic heating, and the 
internal structure is exposed to cryogenic temperature due 
to the hydrogen fuel. Candidate structural design concepts 
are shown in Fig. 10. The paper stresses the importance of 
a highly interactive, multi-disciplinary design/analysis 
procedure. 


Aerodynamic Heating 

The prediction of aerodynamic heating is a problem 
that has challenged thermal analysts since before World 
War II. The problem involves solving the equations that 
describe the conservation of mass, momentum and energy 
for the fluid flow. For altitudes less than 300.000 ft. the 
atmosphere may be assumed to >e a continuum, and the 
set of conservation equations are known as the 
Xavier— Stokes equations. For two-dimensional flow the 
Xavier— Stokes equations governing viscous, compressible 
flow are a set of four non— linear, partial differential 
equations with mixed hyperbolic, parabolic and elliptic 
behavior. A solution of these equations with appropriate 
boundary conditions and constitutive equations for the 
fluid provides distributions for ti e density p. the velocity 
components u and v. the pressure P. and the temperature 
T throughout the fluid. At the fluid 'aerodynamic sutface 
interface, the aerodynamic heat mu i" computed from 
Fourier’s law using the fluid the! nui conductivity and the 
derivative of the temperature normal to the wall. 


flow. These phenomena include: ( 1 ) fluid dynamics effects 
that limit the validity of boundary layer approximations, 
and ( 2 ) high-temperature effects that introduce chemical 
reactions. In addition, at altitudes above 300,000 ft. low 
density effects restrict the application of continuum 
models. 


Classical Boundarv-Laver Predictions 

The concept of the boundary layer was introduced 
by Prandti in 1903 for incompressible flows. The 
extension of the concept to compressible flows began in the 
1930's, and analytical work continued until the mid 1950's. 
Thereafter, most research turned to the development of 
computational methods for solving aerodynamic heating 
problems. 


The boundary layer equations are obtained from the 
Navier-Stokes equations by an order-of-magnitude 
analysis for the thin fluid layer next to the body where 
viscous effects dominate Outside of the boundary— Itiver 
the fluid is assumed inviscicl. and the flow is described by 
the Euler equations. For steady, two— dimensioned, 
laminar flow, the boundary-layer equations are: 
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Prior to the digital computet, t here was no hope for 
solving the complete Xavier— Stokes equations. Instead, 
analysts devoted substantial effort to obtaining 

approximate solutions, and with the help of experimental 
data they developed engineering relations to predict 
aerodynamic heating as well as skin friction. The 
engineering relations typically are based on boundary layer 
or other approximations. The relations are restricted by a 
number of assumptions, and the application of computers 
to the problem of aerodynamic heating has contributed 
significantly to a better understanding of the phenomena 
and improved the reliability of predictions. 

There is substantial literature describing 

aerodynamic heating, and the esoteric details will be left 
for thermal analysts. However, because of the 
interdisciplinary nature of thermal -structural problems for 
high speed vehicles, some fundamental aspects will be 
presented. Early approaches for the prediction of 

aerodynamic heating are described in the 1956 paper by 
Van Driest (ref 24), and in the 1900 text by Robert \Y. 
Truitt (ref 25). A recent discussion appears in the 1989 
text by John D. Anderson, Jr. (ref. 26). 

The preceding section mentioned that aerodynamic 
heating effects become significant at Mach numbers above 
2.5. For Mach numbers up to about five, the flow is 
considered supersonic, and above Mach five, the flow is 
termed hypersonic. The NASP will experience 
aerodynamic heating through the supersonic and 
hypersonic regimes with a peak Mach number of about 25. 
As the Mach number increases above five, new physical 
phenomena become progressively of greater importance 
making hypersonic flow more complex than supersonic 


In th ? above, denotes the pressure at the edge of the 



conductivity is k. and the viscosity 15 p. In general. the>e 
fluid oroperties are temperature dependent. The boundary 
layer equations are nonlinear but exhibit parabolic 
behavior. The perfect gas law is u?od 10 relate predate, 
density and temperature. Fur inn her deiail> including 
boundary conditions see Ref. 2o. 

For some simple but impoitant cases, analytical 
solutions to the boundary layer equations have been 
obtained. The approach involves making changes of 
independent and dependent variables to yield simpler, 
nonlinear, nondimensional equations. These 

translormations lead to the concept of self— similar 
solutions for problems such as a flat plate. Typically, the 
nondimensional equations can be solved numerically for 
various values of the Mach number M. the Reynolds 
number Re, and the Prandti number Pr. 

The solution for the flow over a flat plate obtained 
by Van Driest in 1952 (ref-27 ) demonstrates typical 
featuies. The problem consists of a flat plate with 

specified uniform temperature T , Fig. 11 . The in viscid 

flow outside of the boundary layer has constant velocity 
u x and temperature T . Within the boundary layer the 

viscosity varies with temperature according to 

Sutherland's law (see ref. 26), and the Prandti number Pi 
= /ic /k is assumed constant. Pr = 0.7. The ratio of the 

wall :o free stream temperature is taken as T w /T = 

0.25. Figure ( 11 a) shows typical u velocity profiles, and 
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Fig. ( lib) shows typical temperature profiles- In these 
figures, the Reynolds number, Re x = ^ 3C u 0C x // i xi* 

flat plate boundary layer solution illustrates results 
representative of high speed flows. The velocity profiles 
show that the boundary layer thickness ini leases rapidly 
with increasing Mach number. In fact, Anderson ( ref- 26 ) 
shows that a hypersonic boundaty layer thickness increases 
appproximately as Mach number squared. One 
consequence of a large hypersonic boundary layer thickness 
is that the viscous boundary layer alters the outer in viscid 
flow limiting the accuracy 'of boundary layer predictions. 
This behavior is called viscous interaction. 

The temperature profiles show the fluid behavior 
that is of great importance to the structure. The 
important point is that the peak temperature within the 
boundary layer is higher than the free stream temperature. 
Moreover, this peak temperature increases rapidly with 
increasing Mach number. The high fluid temperature 
within the boundary layer is due 10 viscous dissipation. 
Viscous dissipation is the process where the large kinetic 
energy of the high speed flow is converted to thermal 
energy by the boundary layer shearing stresses. 1 he 
aerodynamic heating rate is proportional, by Fourier's law. 
to the slope of the temperature profile ai the wall (y=o.i. 
Ficure lib shows that the fluid temperatures within the 
boundary layer rapidly min aM j with increasing Mach 
number. Thi> ^ the basic h-mmhi Iui th<* in-vd ol special 
theimal structures tor liinh speni flight. 

Classical, self— similar boundary solution." have also 
been derived for the si agnation region ul a blunt body. An 
important result of this analyse wa> tu show lor 
hypersonic flows that the stagnation heating rate varies 
inversely with the square root o! tne blunt txxly radius. 
Thus to reduce aerodynamic heating the vnnirle nose and 
leading edge regions need to be as blunt as possible 


The Convective Boundary Condit ion 


The heat transfer coefficient is often written in 
non-dimensional form. 


C — 

- p B u p c n 
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or using Eq. (2) 
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H u e c p t 1 aw Tv 


U) 


where C u is the Stanton number Another alternative is 

n 

to define the Stanton number in terms of fluid enthalpies, 
i.e 


C..= 


ws ", nc - K 


) ! 


where H is the adiabatic wall enthalpv. and H is the 
aw * A 

wall enthalpy. In this approach, enthalpv replace' 
temperature as an unknown in the flow analysis Fin a 
c alcrica lly perfect gas (constant specific heats?. Ii = c^T. 

otherwise, dh = c (T) dT. 


Thus, to use a coir. « t live boundary condition it i? 
necessary to know the Stamon number and the adiabatic 
wall temperature. Both. ;;; general, are lottd quantities; 
tha is. they vary with nosition along the wall. The 
classical flat plate bound*: y layer solution will serve as an 
illustration. Extensive studies ol tiie flat plate problem 
fe e ref. 21) show ihar die adiabatic wall temperature cat: 
be ( alculated with good accuracy by 


T = T . 
aw x 



\{ji 


For determining >tructu:al teiupei at tires. aero- 
dynamic heating is olten repieseiucd as a convective 
boundary condition. Boundary layer analysis ol compress- 
ible flow shows that the heat transfer between the 
boundary layei and the wall cat; be expressed aa 

q =h(T -T ) (2) 

*w aw w 

■> 

where q . is the local heating ,aie teg. UTl'/fr - s) at 

the wail, h is the convection coefficient. 1 is the 

aw 

adiabatic wall temperature, and T is the wall 

temperature. The adiabatic wall temperature is the 
temperature the fluid attains for no heat transfer between 
the fluid and wall. Temperature profiles illustrating the 
heat transfer possibilities for high speed flow are shown in 
Fig. 12. The figure shows that for wall temperatures less 

than T , the fluid flow heats the wall: when the 
aw’ 

temperature at the wall is T . there is zero heat transfer; 

and. for wall temperatures greater than T aw , the wall 

heats the fluid. Eq. (2) shows that aerodynamic heating is 
directly proportional to the difference between the 
adiabatic wall temperature and the actual wall 
temperature. The adiabatic wall temperature is always 
higher than the free stream temperature T x and "drives" 

the heat transfer. 


is the ratio ot fluid "pmiii heats * = — and : 

\ 

is the recovery tactor. Analytical and expei imeiital results 
show that the recovery factor can be computed as 


1/2 

laminar flow 

1/3 

turbulent flow 


The boundary layer analysis lor the flat plate flow shows 
that the Stanton number has the tonn 

C H = f (M e , P r , & T w , T e ) !7) 

V Re x 

where f denotes a functional relationship. Eckert in lfloO 
(ref.2S) showed that by using the concept ol a reference 
temperature, a simple formula for Cj.j could be developed 

from the corresponding result for low speed, incompievabh- 
flow. This approach, called the relerence temperature 
method, computes the Stanton number from 


V/if 
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C H = 0.332 ( P r * 


Re.. 


(S) 


where Re x and Pi ate evaluated at a reference 

■ 

temperature T . That is. 


P x 


Re* =-^- 


(9) 




* r> * * 

where P . p and k are evaluated at a reference 

temperature T . The reference temperature is computed 
as 

l’=T [ 1+0.032M J ^0.>> (T /T -I T : 10) 

e [ e v\ ' m ' 

where as before the subscript e denotes quantities at the 
edge of T he boundary layer. Further detail? including 
expressions for skin friction, and equations for turbulent 
How are given in the Eckert papei (ref. 2S). the text by 
Anderson (ref. 26). and other convective heat transfer 
texts. 


Limitations of Classical Boundary Laver Theory 

Classical boundary layer theory provides basic 
insight to understand fundamentals of aerodynamic 
heating. The boundary layer solution in some cases yields 
very practical engineering results that are frequently used. 
But boundary layer theory has limitations that deserve 
mentioning. 


where the flow turns away from the wall into itself 
introducing a local recirculation region next to the wall. 
Flow separation does not occur when the pressure 
decreases in the flow direction. However, when the 
pressure increases in the flow direction (called an adverse 
pressure gradient) flow separation occurs, and boundary 
layer theory does not apply. Flow separation may occur in 
several problems of importance in high speed flight. X-15 
experience provides two such examples: (1) flow over 
protuberances, and (2) shock— boundary layer interactions. 
In both examples an adverse pressure gradient occurs, and 
flow separation develops. A consequence of flow 

separation in both instances is to increase significantly the 
local aerodynamic heating. For these problems, boundary 
layer theory is not applicable, and the phenomena is 
governed by the Navier— Stokes equations. 

Navier-Stokes Solutions 


One of the major advances in computational 
mechanics in the last four decades is the development of 
capability to solve the Navier-Stokes equations 
computationally. The survey article (ref. 30) by Douglas 
L. Dwoyer et al. describes current capability of 
computational fluid dynamics (CFD) for hypersonic 
aircrafi. The Navier-Stokes equations written in 
conservation form are 




tU) 


where L'} is a vector of conservation variables: {Ej and 
{ F } are vectors of the flux comonents in the x and y 
directions. These vectors are 
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Classical self— similar solutions are limited to a few 
problems with simple geometries (e.g. the flat plate J and 
simple boundary conditions (eg constant wall 
temperature). However, in the age of the computer, tins 
limitation is not very serious because the boundary layer 
equations. Eqs. (1), are solved numerically. In fact, since 
the equations are parabolic, inarching methods may be 
used and very efficient computet programs have been 
developed. When combined with in viscid flow programs to 
determine variables at the edge of i he boundary layer, very 
effective analysis procedures have been developed. Some 
of these methods are described in a 1987 paper ( ref .29 ) by 
Fred R. DeJarnette et al. which reviews approximate 
methods for aerodynamic heating analysis. 

There are several flow situations where the 
boundary layer approximation of uncoupling the inviscid 
and viscous analysis does not apply. Two such situations 
are where viscous interaction or flow separation occur. As 
mentioned earlier, viscous interaction refers to flows where 
the viscous boundary layer interacts with the inviscid flow. 
Such a situation may occur for hypersonic flows with 
relatively low Reynolds numbers. Boundary layer 
thicknesses vary inversely with the square root of the 
Reynolds numbers, so low Reynolds number flows have 
thick boundary layers. Flow separation refers to the case 


(F} T = fpv puv 

r 

-io r a 

l \v >’ 


where 8 is the total internal energy. Each of the flux 
vectors contains two vectors of components representing 
inviscid and viscous flux contributions. The conservation 
equations are supplemented by a thermodynamic equation 
of state relating pressure, temperature and density 
(typically the perfect gas law). in the viscous flux 
components, the stresses o . a ami r ate related io the 

velocity gradients assuming Stoke 1 ? hypothesis. The heal 

fluxes q arid q . are related lu lempeiature gradients In 
x y 

Fourier's law. For air, the temperature— dependent 
viscosity is computed from Sutherland's law, and the 
thermal conductivity is computed assuming a constant 
Prandti number of 0.72. Considerable success has been 
achieved in numerical solutions for the two-dimensional 
Navier-Stokes equations with supercomputers. Numerical 
solution schemes include finite difference, finite volume 
and. more recently, finite element methods 
Three-dimensional flows remain a challenge although a 
viscous solution around an X— 24c lifting body (ret. 31) has 
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been obtained. Such solutions, even with the largest 

supercomputers, are not capable ol lesolving 
three-dimensional (low details, particularly aerodynamic 
heating rates. An important pioblem for which 
Navier-Stokes solutions are providing valuable insight is 
t wo _di m ensional shock interaction phenomena for the 
NASP. 

In 196" NASA conducted a series of X— 15 flights 
with a dummy hypersonic ramjet engine mounted on a 
pvlon under the rear of the fuselage. On the third night 
with the dummy engine on October 3_, 1967 the X-15 
reached a maximum Mach number of 6.7 at an altitude of 
99.000 feet. During the flight severe structural damage 
was experienced due to complex shock impingement and 
interference effects on local aerodynamic heating (ref. 32). 
Considerable heating— induced damage occuned on the 
engine pylon indicating that local temperatures exceeded 

the Inconel X melting temjierature of 2600 F. Since then, 
shock interference heating has been recognized as a critical 
problem for high speed vehicles because extreme pressure 
and heat transfer rates tan omit in highly localized 
regions where the interference pattern impinges on the 
surface. 

Shock interference heating is an important 
consideration in the design of the cowl leading edge of the 
engine structure of the NASP The problem has been 
strong motivation for recent studie> ol shock interference 
heating on leadins edges. Figure 1-3 shows the overall flow 
field for the NASP. The interaction between the aircraft 
bow shock and the cowl bow shock can cause se\eie local 
aerot hernial loads on the cowl leading edge Shock 
interference heating on leading edges has been studied 
experimentally by Wieting and Holden (lej. 33) and 
tompututionailv by Klopiei and \ w ret. 34-3o.i as well as 
Thareja et al/ (ref. 36). Dethaumplui ei al (ret. 3<) 
studied the flow-thermal-strucuual behavior of leading 
ed*es. Figure 14 shows the problem statement for the 
Navier-Stokes solution, an adaptive finite element mesh 
and a comparison of experimental and predicted heating 
rate distributions. Adaptive, unstructured meshes, a 
recent advancement in finite element and finite volume 
flow analysis, provide high resolution ol flow details such 
shocks and boundary layers "1 he comparative :>uitace 
heating rates shown in Fig. 14c >how very good 
agreement, but the paper points out there was some 
difference in prediction of the undisturbed heating rate q Q 

that normalizes the curves. Note that the results show the 
interference heating rate is almost 15 times the 
undisturbed level. Peak interference heating rates for 
NASP flight conditions and geometries can be as high as 

70 000 Btu/ft 2 - s. To put this value in perspective, a 
typical 1 kW portable hair dryer produces a heating rate ot 

about 30 Btu/ft 2 - s. Such high local heating presents a 
formidable challenge to the structural designer. 

A discussion of aerothermodynamics of 

transatmospheric vehicles from the design perspective 
appears in the papers of Tauber et al. (refs. 33—39). 

Aerothermal Load Effects on Flight StruCWS S 

In his classic 1956 paper (ref. 6), Professor R. L. 
Bisplinghoff identified the basic structural and aeroelastic 
considerations for high sj^eed flight. From a structural 
perspective the considerations include, (i) deterioration of 
mechanical properties at elevated temperatures. (-) 

thermal stresses introduced by temperature gradients. (3) 
modification of stiffness and vibration properties, and (4) 


aeroelastic instabilities. These considerations remain 
important todav. The general problem consists of 
determining the behavior of the flight structure under 
aerothermal loads, primarily aerodynamic pressures and 
heating. The behavior of the flight structure can be 
classified according to whether the response is quasi-static 
or dynamic. In a quasi-static response inertia forces are 
negligble, and the structure responds to the aerothermal 
loads slowly. In a dynamic response inertia forces have a 
significant role, and the structure responds with vibratory 
motions. The determination of the quasi-static response 
of structures has traditionally been called theimal stiess 
analysis. The determination of the structural dynamic 
response considering the interaction of the deforming 
structure and the aerodynamic pressure in the presence of 
aerodynamic heating is called aerothermoelasticity. An 
area of recent concern is the response of structures to 
acoustic loads. The determination ot the dynamic response 
of a^ structure to random fluctuating pressures in the 
presence of aerodynamic heating may be called 
aerot h er moacoust i cs . 

Quasi -Stat ic Interactions 

To determine the quasi— static response, a logical 
appiodch is to sepai ate the aerot hermal-smictuial 
problem into distinct uncoupled problems by assuming 
weak coupling between the external aerodynamic flow and 
the $ ruttural response. This assumption is permissible 
when structural deformations are too small to alter the 
external flow. During X — i "» flights irels. 9-1U) seveial 
unexpected thermal problem* wen* encountered due to 
intense local aerudvnamic imuiing because ol quasi-static 
interactions. Windshield damage occurred when thermal 
buck! inti of the retainer liamt- uimed niLeu>e local heating 
in th? class. These pioblems were identified and solved 
during flight tests. However, accoiding to Ret. 10. the 
realK important lesson here is that what are ini not and 
unimportant features ot a subsonic ui supersonic anctaft 
must be dealt with as prime design problems in a 
hypersonic airplane". 

A more recent example of a quasi-static 
flow /structural interaction is thermal protection sv stems 
tested in the Mach 7 eight-foot High Temperature Tunnel 
at the NASA Langley Research Center (ref. 40). The tests 
show that panels "bowed— up" into the flow to pioduce 
heating rates that are up to 15 times greater than flat 
plate predictions. Thornton and Dechaumphai (rei. 41) 
used a finite element approach to study coupled flowv 
thermal and structural behavior of aerodynaniicaily heated 
panels. Some typical results lor concave and convex 
deformations are shown in Fig. 15. Only very modest 
deformations occuned. but flow teatures were altered 
significantly. For the Mach 6.6 conditions studied, panel 
deformations introduce shocks, expansions and 
recirculation regions in the flow. The effect of convex 
panel deformation is to increase local heating rates on 
windward surfaces and decrease them on leeward surfaces. 

In the overview paper (ref. 23) describing the design 
of an airframe structure for NASP, the interaction between 
the external flow and the thermoelastic deflections ol a 
movable wing is described. Figure 16 shows the results of 
the aerot hermoelas tic load analysis. The wing deflected 
sha]>e is shown for several different iterations in the 
analysis cycle. The paper notes the consideration of the 
interaction resulted in the prediction of lower deflections 
and stresses which translated into reduced structural 
weight. * 

The recent analyses of quasi— static flow /structural 
interactions illustrate some of the important ellecta ot this 
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behavior on flight structures. However, because of their 
multi-disciplinary nature, the analyses are difficult and 
expensive. In addition, the analyses are unvalidated by 
experimental data. 

Dynamic Interactions 

In a 1963 paper (ref. 42) I. E. Garrick of NASA 
Langley surveyed developments in aerothermoelasticity. 
The classical aeroelastic triangle representing interactions 
between the fields of aerodynamics, elasticity and inertia 
was extended to include thermal effects. The new figure, 
the aerothermoelastic tetrahedron, illustrates the 
interdisciplinary aspects of aerothermoelasticity. The 
paper also discusses aeroelastic consideration of the X— 15, 
effects of transient heating on vibration frequencies, and 
panel flutter. Early in the X — 15 flight program t lie pilot 
reported a rumbling noise at high dynamic pressures. The 
turned out to be panel flutter ui large aieas ol the skin ol 
the side failings and tail. This problem was strong 
motivation for NASA studies in the 1960's and 1970's of 
panel flutter at elevated temperatures. A recent 
assessment (ref. 43) of flutter model testing relating to 
NASP provides an excellent summary ol flutter literature 
for supersonic and hypersonic flight. 

The response of aerodynamically heated structures 
to acoustic loads is an important consideration because of 
acoustic latisue. Rehuence 44 notes that operational 
experience for a variety of airciaft has demons! rated that 
intense acoustic pressures can cause fatigue failures of 
lightweight structures Preliminary estimates of acoustic 
loadings for the NASP indicate pressuie levels that are 
well into the range where acoustic fatigue failures have 
occurred in the past. Thermal effects are important 
because: ( 1 } there is little acoustic data for high 

temperature structures and materials. (2) thermal 
prestrain and buckling can aflec strain levels significantly, 
and (3) temperature can change material properties and 
the latiuue lifetime as expies?ed in the S-N diagram. 
Recent efforts lo develop toinput.tuoinil method* lor 
aerothennoacoustics aie destiibed in the dissertation by 
Locke (rel. 45) and the paper bv Locke and Mei (ref. 46). 

Design of Thermal Structures 

The design of thermal structures is a complex 
process that involves consideration of the flight vehicle 
trajectory, aerothermal luads. thermal structural concepts 
and materials. A near optimum design involves tradeoffs 
among these and other factors. The design of thermal 
structures is too complex to be discussed in detail in this 
paper, but fundamentals will be cited. 

Fli ght Regime 

To determine aerothermal loads and heat transfer 
to the vehicle, the flight trajectory must be determined. 
The trajectory is determined by mission requirements. 
The basic equations of flight mechanics, hypersonic 
aerodynamics and re-entry heating appear in the text by 
Wilbur L. Hankev (ref. 47). Recent design studies for 
hyper velocity vehicles (ref. 4$, and for NASP (ref. 23) 
describe relationships between mission requirements and 
thermal-structural design concepts. 

A parameter used in the design of thermal 
structures is radiation equilibrium temperature. The 
radiation equilibrium temperature is the upper level that 
the surface of a structure can reach. An energy balance at 
the surface of a structure at radiation equilibrium states 
that the aerodynamic heat flux given by Eq. (2) is equal to 


the heat flux emitted by radiation, i.e . 

<1. - « IT..-V - « V 1131 

where T f is the radiation equilibrium temperature. < is the 

surface emissivity, and a is the Stefan-Boltzman constant. 
Eq. (13) states that all of the incident aerodynamic heat 
flux is emitted by radiation; none of the incident flux is 
conducted into the structure. Actual surface temperatures 
ma> be lower than radiation equilibrium temperatures 
because of conduction heat transfer into the structure. 

Two key factors of the flight regime for thermal 
structural design are radiation equilibrium temperatures 
and exposure times. Radiation equilibrium temperatures 
for the X— 15, the Space Shuttle, and the X-30 can be 
compared in Fi|s. 3, S and 9. respectively. The flight 
regime for NASP vis-a— vis other flight vehicles is shown 
in Fig. 17. The important point is the design of hypersonic- 
cruise vehicles like NASP differ from previous designs 
because of much longer flight durations. A third 
important factor of hypersonic flight is use of liquid 
hydrogen fuel. Among other considerations, cryogenic 
hydrogen introduces significant smicmral temperature 
giadioius. Both maximum tem;x*i at 1 1 res and gradient.' 
influence structural design concept rind material selection. 

Thermal Protection Systems 

Hypersonic flight vehicles require thermal 
protection systems to withstand sustained aerothermal 
ioac?. The thermal protection system, the supporting 
diriiame and the engine structure are examples of thermal 
structures. The design of <\ thermal protection system is 
bated upon the principle that the energy transmitted by 
I he hut boundary la>er flow must be absorbed or rejected 
by the thermal piotection system. Figure IT identifies 
theimal protection systems used by hypersonic flight 
veh ties. For relatively short missions the Mercury and 
Apollo spacecrafts absorbed the thermal energy through 
ablative heat shields. The Space Shuttle absorbs re-entry 
heading and thermally insulates the airframe by a very 
effective, but fragile, tile system. For relatively short 
missions, the X— 15 absorbed the aerodynamic heating by 
using skin and airframe of high temperature metallic 
materials. For long duration but lower temperature flights 
the VF— 12 (or SR — 7 1 > u>ed a high— temperature titanium 
st met lire with high emissivity >ui lace coating at ladiatiou 
equilibrium to reject the aerodynamic heating. For 
sustained, very high temperature flights a combination of 
several thermal protection concepts is likely to be used. 
Figure 10 shows candidate thermal piotection systems for 
NASP which include hot structures, insulated structures 
and convectivelv cooled structures. Convectively cooled 
airframe structures are relatively new with no previous 
flight experience. Design and experimental studies (ref. 
IS) conducted over the last twenty years indicate that 
these systems will be effective in absorbing thermaJ energy 
in regions of intense local heating, e.g. shock interference 
heating on engine structures. 

Ma erials 

The evolution of thermal structures for high speed 
flight vehicles has placed ever-increasing demands on 
material performance. Structural designs typically require 
high-stiflness, thin-gauge materials that can be fabricated 
into complex, built-up structures. Thermal structural 
designs typically require high strength, low density 
mai erials that retain i heir desirable properties at elevated 
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temperatures. There is a substantial research effort 
currently underway to develop new aerospace materials to 
meet these challenges. References 49—56 are recent papers 
describing progress in material development for thermal 
structures. 

The performance of the NASP will depend on the 
development of new, lightweight materials that can 
perform at temperatures higher than today's materials can 
stand. Today’s high-temperature materials such as 
nickel-based alloys cannot be used because structural 
weight limits mandate lighter-weight materials. The 
structural weight limits require the use of thin gauge, low 
density alloys and composites. Figure IS compares specific 
strength and stiffness of several advanced materials. 

Below 2000 °F the structural designer has several 

alternatives, but above 2000 °F choices are limited. The 
NASP materials and structure augmentation project (ref. 
22] lias the specific goal of advancing the 
readiness-for-use data of several advanced materials that 
are needed for the vehicle. The contractor and materials 
under development in this program ate listed in Table 1. 


Table 1 

NASP Materials and Structures Augmentation Program 


Contractor 

Material 

General 

Dynamics 

Refractory composites 

carbon-carbon composites 
ceramic composites 

McDonald 

Douglas 

Titanium metal matrix composites 

Pratt 

and Whitney 

High creep strength materials 

Titanium aluminide alloys 
Titanium aluminide composites 

Rockwell 
North American 

Titanium aluminide alloys 

Rockwell 
Rocket dyne 

High thermal conductivity composites 
Coppei matrix composites 
Beryllium alloys 


Heat Transfe r and Thermal Stresses 

The design of thermal structures requires the 
determination of temperatures, displacements, stresses and 
strains throughout the structure. A structural heat 
transfer analysis is needed to determine maximum 
operating temperatures that guide material selection. 
Temperature distributions are also needed to compute the 
" thermal— loads" for a structural thermal— stress analysis. 
The heat transfer and structural analyses rest on the 
conservation equations of continuum mechanics, 
constitutive models of material behavior, and 
computational methods implemented on modern 
computers. 


Conservat ion Equations 

For continuum formulations of solid mechanics, the 
conservation equations that must be considered are 
conservation of linear momentum, conservation of angular 
momentum and conservation of energy. The conservation 
of linear momentum produces three equations of motion 



+Bj=/> 


d 2 u. 

d T 


( 14 ) 


where a - are components of the stress tensor, IT denote 

bodv orce components per unit volume, p is the material 
density and u- are the displacement components. 


Conservation of angular momentum shows that the stress 
tensor is symmetric. Conservation of energy considers the 
work done by the stresses, thermal energy transported 
across surfaces by conduction, thermal and mechanical 
energies stored within the liuitei ial. *ind kinetic enetg\ due 
to material motion. The conservation ol energy equation 
for a (onrinunm model of a delormable ImkIv is 


d q 
d xi 


d (■; 


M ur 0 u n 

xT^ijFt+'FT-Q 
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where c\ { denotes components of heat flux, are 

components of the strain tensor, u is the internal energy 
per unit mass, and Q is the mte of internal heat generation 
|XM' unit volume. The internal energy ol the solid depends 
on 1 1 u* Stiains ami icmixu me. that i?. u = u i T). 

The strain component > at a point are related to the 
di'ph cement components b> 


, u- , u 
d i , d j 

+0 

"k d V; 

0 X d X; 


x . d x • j 


Ofter displacement gradients are assumed small, and the 
last term in Eq. (16i that involve? pioducts and powers ot 
displacement gradients is neglected in comparison to the 
first two terms. Tin* result is the linear 
rail .—displacement relations 
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Equations (16) or (IT) show that the strain tensor is 
symmetric. Thermal-structural problems are olten 
formulated in terms of the linear strain— displacement 
relations. However, under severe conditions, structures 
may experience large deformations, and the nonlinear 
strain-displacement relations may be required. In these 
circumstances, the definitions of the stress components 
have to be— interpreted relative to the undeformed and 
defunned configurations of the structure. The 
conservation equations, in the context of thermal stresses, 
are derived in the classic text by Boley and Weiner (ref. 
57), and in the recent series of volumes on thermal stresses 
edited by Hetnarski (refs. 58-59). 

Conservation of energy. Eq. (15), states that there 
is a relationship between stresses, strains and temperature 
in a deformable body. In one interpretation, conservation 
of energy indicates that variation? of stresses and strains 
within the solid alter the heat flow and thermal energy. 
The equation expresses a conversion of mechanical to 
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thermal energy. The relationship between stresses, strains 
and temperature in the energy equation is known as 
thermal-mechanical coupling. The conversion of 

mechanical energy to thermal energy according to 
conservation of energy is a well-known phenomena which 
has been studied extensively for elastic behavior (ref. 60). 
Studies of coupled thei moelasiiuty have shown that for 
metallic materials within the elastic range, the conversion 
of mechanical energy to thermal energy may be neglected 
for aerospace applications. More recent studies (refs. 
61—63 ) of thermal— mechanical coupling for non-linear, 
inelastic behavior suggest that coupling may be important 
under some circumstances. However, in analyses of flight 
structures, thermal-mechanical coupling is usually 
neglected. One argument for this assumption is that the 
external energy supplied to the structure by aerodynamic 
heating is so large that the thermal energy converted from 
mechanical energy is negligble in comparison. This means 
that the energy equation can be simplified by assuming the 
solid is undeformable. that is <^=0 which is the usual 

approach in heat transfer texts. Tndei this circumstance 
the internal energy is regarded as a function of 
temperature alone, and it is customary to take 
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where t(T) is the mat‘Mirtl'> >pr(ilii heal which i> 
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which is the equation customarily found in heat transfer 
texts. The heat flux component? aie normally related to 
temperature gradients by Fourier's law. For an 
anisotropic material Fourier'? Ian ?witev 


q i --k ijT (19) 

1 j 

where k- j are components of a thermal conductivity tensor. 

In general, the material thermal conductivities are 
temperature dependent. 


The heat conduction equation is solved subject to 
an initial condition and boundary conditions on all 
portions of the surface. The initial condition specifies the 
temperature distribution at time zero. The boundary 
conditions may consist of specified surface temperature, 
specified heat flow, convective heat exchange and radiation 
heat exchange. These may be written as 


T S -T 1^ X 1, X 2, X 3,^ 

on Sj 


Qi n--q s 

On S ry 

(21) 

Vr h (T s — T e ) 

on 


Qj n i =«Tj-Qq r 

on S 4 



where n . denote components of a unit outward normal, and 
S| (i—1.4) denotes portions of the surface. The specified 
surface temperature is and the specified surface heat 
flux is q g (positive into the surface). In the convective 

boundary condition the convective exchange temperature 
T e is T (see Eq. (2)) for aerodynamic heating. In the 

radiation boundary condition q is the surface absorptivity, 
and q is the incident radiation heat flux. For structural 

heat transfer with surfaces at significantly different, 
elevated temperatures, radiation exchanges between 
surfaces can occur. The determination of radiation 

exchanges between surfaces is complicated because:! 1) 
radiation emitted by a typical surface depends on its 
surface temperature which is unknown, and (2) the 
geometrical relationship between each surface must be 
considered. Often, the radiation problem is handled by 
discri t izing the radiation boundary into N discrete surfaces 
which are assumed isothermal. If a radiation heat flux on 
the it h surface is called H-, a set of N simultaneous 

equations may be developed to determine Hj in terms of 
the temperatures T«. In matrix form these equations may 
be written as 

[[I] - [F] (*-<)] (H} = [F] {fffT 4 } (22) 


Structural Heat Transfer 

For conduction heat iran.dn in .irtospace 
structures, the classical hc.ii coinluciion equation i> used. 
Substituting Eq. (19) into Eq. tl>) yields 



where the components of the matrix [F] are the 

viewfactors F^. [I] denotes the identity matrix, and is 

the Faction of the radiation energy leaving surface i that 
arrives at surface j. The determination of viewfactors for 
complex three-dimensional structures is a formidable 
computational task. 


which is a parabolic partial differential equation. This 
means that thermal disturbances propagate at infinite 
speeds through the body. To address this anomoly, some 
authors (refs. 60, 64-65) have used a modified form of 
Fourier's law to derive a hyperbolic energy equation. In 
the hyperbolic energy equation, thermal disturbances 
propagate at a very high, finite wave speed which is called 
"second sound". The application of the hyperbolic energy 
equation for solids is controversial because the phenomena 
of finite wave speeds of thermal disturbances has never 
been demonstrated for structural materials although it has 
been for gases. The parabolic energy equation is 
customarily used for structural heat transfer. 


An excellent example of heat transfer for a complex 
structure with aerodynamic heating is the space shuttle 
wing. William L. Ko and co-authors have documented 
several studies (e.g. refs. 66-68) of heat transfer analyses 
with comparisons to flight-measured temperatures. The 
wing geometry, a three-dimensional model of a wing 
segment, and comparative histories of structural 
temperatures axe shown in Fig. 19. The calculated 
aluminum structural temperatures agree reasonably well 
with the flight data from re-entry to touchdown. The 
effects of internal radiation and internal convection were 
found to be significant. 
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For conveciivelv tooled stiuctures the heal tiansfei 
in coolant passages must be considered. The dominant 
nn>de of heat transfer in the coolant How is forced 
convection. An engineering model ot flow in the coolant 
passage is typically used. The engineering formulation is 
based on assumptions that produce a one dimensional 
energy equation with the bulk temperature T^x,t) of the 

coolant as the fundamental unknown. The coolant energy 
equation takes the form 

- l*t k f A f§^ +mc f fr - h * T w- T f>^f c fw =0 (23) 

where the subscript f denotes fluid quantities. In the 
above Aj- is the cross-sectional area of the coolant passage. 

rh is the coolant mass flow rate, h is a convective 
coefficient describing heat exchange between the wall ol 
the coolant passage and the coolant, and p is the coolant 
passage perimeter. Structural tempera lures and coolant 
temperatures are determined by solving the energy 
equations. Eqs. (20) and l21). simultaneously Fuithet 
detail.- of the formulation and numerical example- are 
described by Thornton and \V leiiug m dels. lit)-. i. 


The equilibrium equations are solved subject to boundary 
conditions specifying, either the displacement components 
or surface tractions on all external surfaces ot the 
strut-inn* The initial t onditioiis include specifying values 
of the structure's displacements, and if the structure has 
lieen .abject to previous loads, initial values of stiesse? anu 
strains may be required 

For linear, elastic behavior the stress components 
are elated to the strain components by generalized 
Hooke's law. For a homogeneous, isotropic material this 
constitutive relation may be written as 

< 7 . . =: A + 2 G { 3 A+2G ) ^. a(T-To) (25) 

i) ij kk ij *J 

whore A . is the Kionet ker deb a; X and G are the Lame 
constants; a is the coefficient of thermal expansion, and l u 

is the reference temperature for zeio thermal stress. The 
Lame’ constants are related to more familiar engineering 
run? ants by 


Over the la.-t four decades there have been 
significant advances m computational methods tor 
structural heat transfer. The use ul computers and 
computer graphics has made analyses of complex theimal 
structures a routine step in the design process. The 
proceedings iref. r2j ot a conference on computational 
aspect? of heat transfer in sumctuies held at NASA 
Lanaley in 10S 1 describes capability to compute 

temperatures ut flight vehicles. Cumptitei haidwaie and 
software have made and continue to make significant 
advances. So It ware is based on the linite 
diffeiente/liimped— parameter method or the finite element 
metliud Program? ba?ed on tne toimei meihuti include 
TR AS VS. MITAS and SINDA. TRASYS has been used 
extensively for L. S spacecraft radiatiuii heat uanster 
analysis since 1972. There are several widely-available 
finite element programs with heat transfer capability 
including ABAQLS. ANSAS. EAL. MARC. 
MSC/NASTRAN and PAT R AN. Current computer 
hardware trends include iinpiuvmg cumputatioual .'peed? 
by vector and/or parallel processing 


Thermal Stresses 

Thermal stresses in a structure are determined by 
solving the equations of motion, Eq. (14) and the 
strain— displacement relations. Eq. (16) simultaneously 
with constitutive equations relating the stresses and 
strains. One of the basic assumptions in thermal stress 
analvsis of flight structures is that the inertia forces in Eq. 
(14) can be neglected. In most practical applications the 
thermal response of a structure is relatively long in 
duration compared to characteristic times of the structural 
response. Under these circumstances, the structure 
responds to a time— varying thermal load in a quasi— static 
manner proceeding through a succession of equilibrium 
positions without oscillations. Thermally induced 
oscillations occur only when the thermal response time and 
the structural response time are about the same. Suddenly 
applied thermal loads of high intensity to thin beams, 
plates or shells can induce oscillations (see ref. 57), but 
typical flight structures respond tu aerodynamic heating 
quasi-statically. Neglecting inertia forces means that in a 
quasi-static thermal stress analysis the equations ol 
motion reduce to the equilibrium equations. 


e, “2< ltn 

whe e L h the modulus ul elasticity, and r is Pui?-pn * 
uitu> for small leinpeialuu change? \ lie elastic ptopertn- 
are ;onstant. but for large temperature changes the elasi u 
properties are temperature-dependent. Thus il the 
temperature varies throughout the structure, the 
properties vary from point to point. For small stresses 
strains and/or temperature changes, the behavior ul 
structural members is elastic, and the solution ot 
boimcUy/initial value, problems makes up the Held ot 
ihemioeiasiicity (ref. 57-bO). 

One of the most signiticant developments ol the las ^ 
foui decades is the finite element method. M itii linite 
elements, virtually any complex structure can be modeled 
and analyzed with a high degree of accuracy. LiteialU 
dozens of finite element books exist and a vast literature 
comprising thousands of papers describe finite element 
me i hodology. Although the method originated for aircratt 
structural analysis, finite elements enjoy success in related 
thermal— structural disciplines including heat transtei and 
compressible flow analysis. Previous sections have 

iiRMiiuiicd these application ami given ?elecied u-teieuce' 
Finite element thermal stress analysis capability e\i>i- in 
varietv of commercial software including the uku- 
mentioned in the preceding section on structural neat 
transfer. 

Under high rates of loading with the material at 
elevated temperatures, flight structures will experience 
inelastic behavior that includes rate-dependent plastic 
(viscoplastic) deformations. Most metals exhibit 
viscoplastic behavior at temperatures above 40 7c oi then 
melting temperature. One of the significant advances in 
the last twenty years is the development of constitutive 
models that include both plasticity and creep in a single 
set of equations called unified constitutive equation? 
Several investigators have developed unified constitutive 
models. The models are based on microphysical benavioi 
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of materials. are guided by phenomenological 
considerations, and employ concepts of continuum 
mechanics. Typically, the equations use the concept of 
internal state * variables to represent the evolution of 
material behavior. The equations also involve a number of 
material parameters (some of which ate temperature 
dependent) that must be determined experimentally. 
Reference 73 contains articles describing several 
constitutive models. 

One of the most well known of the constitutive 
models was introduced by S. R. Bodner in 1968. A survey 
article in Ref. 73 by Bodner describes the evolution of his 
unified constitutive model. The approach has become 
known as the Bodnei-Parton constitutive model. To 
consider the time dependent character of t he behavior, the 
equations are written in rate lotm The loial strain rate 
tensor is separated into elastic and inelastic components. 

E 1 

i \ 


where a dot denotes differentarion with lespect to time. 
The elastic component indicated by the superscript E is 
related to stresses by the time derivative of Hooke's law. 
Ec|. (23). The inelastic component is derived beginning 
with the Prandtl-Reuss flow rule of plasticity. The 
inelastic strain component is represented by Boclner and 
Parton as 


wide temperature range. Material property data for the 
Bodnei-Parton model is available for Ti, Cu, Al. Rene' 95, 
IN- 100, Inconel 7 IS. Hasteiloy-X, Bl900+Hf and 
Mar-M247. 

Unified constitutive models have been implemented 
into finite element analysis by a number of researchers and 
in t ie commercially available MARC program. The paper 
(ref 75) by Thornton et al. contains references to these 
finite element viscoplastic analyses and describes 
viscoplastic analysis of hypersonic structures subjected to 
severe aerodynamic heating. A thermoviscoplastic analysis 
was performed for a convectively cooled segment of the 
scram jet engine structure for the NASP. The quasi-static 
thermal-structural analysis was performed with a finite 
element mesh (Fig. 20a) for a nickel super allov 
repiesented by the Bodner-Parton constitutive model. 
Figure 20b presents of the temperature history on the 
aerodynamic surface indicating the very high temperature 
the material experienced. Figure 20c presents the history 
of the normal stress component at the aerodynamic 
surface. Note the material yields early in the response. 
After the temperature drops there is a rapid decay ol the 
high compressive stresses, and the analysis predicts 
residual tensile Stress in the plastically deformed suitac** 
Recent thermal— struci ural analyses of hydrogen— <. ooled 
leading edge concepts for hypersonic flight vehicle are 
lies :ri bed in refs. 76-77. 

T henna [-Structural Testing 


1 s. 
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where S is the deviatoric stmss tensor, and is the 
U . - 

-ecund deviatoric sues.' invariant The coeiiiciciu is a 

material parameter, n is a temperature dependent material 
parameter, and Z is an unknown load history dependent 
parameter called the internal state variable. I he evolution 
equal ion lot the slate variable lor inotropic naidening is 


\i- Z.,l'l 

2=mjZ,— ZAV.-A.Z. 1-^-p (29) 

11 1 Up! j 

with the initial condition Z(o)=Z o In this equation, m^. 
7.y and Z 0 are specified parameters: .-Q and ij are 

temperature dependent parameters: aud \Vj is the inelastic 
work rate, 

W I =<7 ij‘ij (30) 

In contrast to classical elasticity and plasticity, the 
Bodnei-Parton constitutive model assumes that lor all 
load levels there is inelastic strain. In addition, the model 
does not employ a yield criterion, in applications, the 
inelastic strain component is small in comparison to the 
elastic component at low load levels and becomes 
significant only when inelastic phenomena become 
prominent. A NASA-Lewis S{)onsorecl research program 
(HOST) conducted by the Southwest Research Instil utie 
recently concluded a four year effort (ref. 74) to lutther 
develop unified constitutive models for isotropic materials 
and to demonstrate their usefulness for analysis of high 
temperature gas turbine engines. One result of this study 
is material property data Tor nickel-based alloys over a 


Thermal-structural testing remains an important 
step in the development of .structures for high speed flight. 
Test facilities exist at several NASA and Air Force 
ins .dilations as well as in private industry, but the 
principal l'. S. government facilities are at the Air Force 
Flight Dynamics Laboratory (ret. 7 S ) and the NASA 
Di deu Flight Research Facility (ref. 79). Two recent 
Mudies related to the NASP are good sources of 
information about past test programs and current 
capabilities. In Ref. SO. Inger P. Friedman et al. provide 
an assessment of thermoplastic analysis ancl testing 
applicable to NASP: in Ref. 81. H. A. Hanson and J J. 
Ca>e\ describe a study to determine high temperature 

( 1 COO — 3000 °F) capability for testing full-scale aerospace 
vehicle structures. Over the years NASA Dryden has been 
in\ol\ed with extensive flight and laboratory test 
programs for the X— 15. the N T— 12. the Spate Shuttle and 
hv :>ei>onic structural components 

An excellent review of hot structures test and 
annlvsis technology is the proceedings (ref. 82) ol a 
wuikshop held at NASA Dryden in November. 1988 In 
summary V. Michael DeAngeiis, the workshop coordinator, 
made the following observations: (1) hot structures testing 
is expensive requiring sophisticated computer control 
systems, large amounts of instrumentation, and high power 
requirements, (2) it is time consuming because of the need 
for extensive instrumentation checkout and calibration, 
and (3) test procedures are in an early stage ol 
development for new, high temperature structures. He 
noted that correlation of test data with analysis is 
becoming more difficult as Mach number increases 
because: (1) structures are becoming more complex with 
new materials and active cooling, (2) computational 
complexity is increasing with the need for finer models to 
capture thermal gradients and high local stresses. 13 j rest 
requirements are increasing, e.g. more instrumentation is 
needed. (4) measurement a capability decreases with 
increasing temperature; e.g. accurate strain measurements 
are particularly critical, and (5) test capability decrease:- 
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significantly at high temperatures. 

A conclusion that may be drawn from reviewing 
this recent literature is that there is a significant, real need 
to develop new high temperature test technology. This 
technology is needed to assess the performance of new 
materials and design concepts as well as to validate the 
analysis tools required to develop thermal structures for 
high speed flight.. 

Concluding Remarks 


This paper surveys progress m thermal structures 
from the early days ol supersonic flight to the cuneni 
research and development for the National Aerospace 
Plane Fundamental concepts of aerodynamic heating, 
aerothermal load effects on flight structures, design of 
thermal structures as well as heat transfer and t hernial 
stress analysis and testing are reviewed. Major advances 
in technology have occurred that provide the foundations 
for the design of thermal strucuues lor flight vehicles in 
the twenty— first century. Much progress has been 
accomplished, yet there are a number of research needs 
i kit must be addressed: 


The X— 15 was the first, 
manned vehicle capable 
missions at Mach 5 at alu 
higher. The Iasi X— 15 fl 
X— 15 Higlus were enomu 
maiiv significant comributi 1 
of hypersonic flight nn Inch 
structures Although t 
( unu ibuted >iimili< ant iy 
program. T l n • i * - ' ei 

• experiment al. ii> poi Mjnn tl 
i he world leader in high 
must not laltei in elloit? to 


and to date, the only 
of living atmospheric 
nudes of 100.000 feet or 
“lit was in 1 *)6S- The 
ni^ly successful making 
mis to the understanding 
m» the design of thermal 
ha Space 'Shuttle has 
1 1 1 i In- nation s space 
v *ai<)H“ mvd tin an 
iiihi \ » *1 1 k le. f o icmain 
npiH'd t Sight . the nation 
develop t lie N ASP 


Significant advances liave occurred in 
computational fluid dynamic? and computet 
hardware that permit high (piality solutions lor the 
Navier— Stokes et|uations. Accurate prediction of 
aerodynamic heating lor two— dimensional flows is 
possible, but the prediction ol aerodynamic heating 
in three-dimensional flows lemains a challenge. 

The importance ol interactions between high speed 
Hows and hot. deforming strucuues has been 
recognized. but analyses ot coupled 
flow— thermal-structural interactions is in an early 
stage of development. Computational studies ot 
interaction effects should continue. There is also a 
need for experimental data to support 
computational studies of interaction eflects. 

Substantial progress lias been made in the design 
and development of con vet lively cooled structures 
for high speed flight, but there is a clear need lor a 
flight test program to validate the designs under 
realistic conditions. 

The need for new, lightweight materials with 
well— understood behavior for structural 

applications above 2000 °F is critical. High 
material costs currently limit experimental studies 
particularly in university research programs. 
Substantial basic research studies of new materials 
such as metaJ matrix composites are needed. 


Thermal-structural analysis capability with the 
finite element method has reached an advanced 
stage of development. New developments in 
constitutive modeling permit analysis of highly 
nonlinear material behavior. For complex night 
structures, analysis capability has exceeded the 
available experimental data base. There is a need 
for high temperature experimental data to validate 
analysis capability. 


New high temperature test technology' is needed to 
support experimental studies of new material?* and 
design concepts at elevated temperatures. 
Development of methods for accurate measurement 
of strains at elevated temperatures is a high 
priority. 


Manv of these recommendations should and are 
being pursued to support the development of ihe NASP 
In mid it ion. the nation must maintain ^ h^ad-basei. 
fundamental research program. Utchard R Helcimle.' 
looking back in 1962 on NAT A thermal-suumiia! 
research from 194$ to 1956 made the argument tor basic 
research very well: "A healthy research program must 
provide freedom to explore new ideas that have no 
obvious applications at the time. These ideas may 
venerate the technology that makes impoitant. 


* 
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